##Gerando dados
set.seed(42)
area <- rnbinom(10,mu=500,size=0.5)
riqueza <- rpois(10,lambda=(area^0.38))

##Tutorial
area <- c(300, 350, 961, 295, 332, 47,  122, 11, 53, 2749)
riqueza <- c(1, 7, 20, 7, 8, 4, 8, 3, 5, 23)
area
riqueza
summary(area)
summary(riqueza)
mean(x=area)
varea <- var(area)
varea
sqrt(varea)
sd(x=area)
mean(riqueza)
var(riqueza)
sd(riqueza)
plot(x=area, y=riqueza, xlab="Area (ha)", ylab="Nmero de Espcies")
modelo1 <- lm(riqueza~area)
summary(modelo1)
previsto <- fitted(modelo1)
riqueza - previsto
residuals(modelo1)
par(mfrow=c(2,2))
plot(modelo1)
par(mfrow=c(1,1))
plot(x=area, y=riqueza, xlab="Area (ha)", ylab="Nmero de Espcies")
abline(modelo1)
plot(x=area, y=riqueza, xlab="Log Area (ha)", ylab="Log Nmero de Espcies", log="xy")
modelo2 <- lm(log(riqueza,base=10)~log(area,base=10))
par(mfrow=c(2,2))
plot(modelo2)
par(mfrow=c(1,1))
plot(area, riqueza, xlab="Log Area (ha)", ylab="Log Nmero de Espcies", log="xy")
abline(modelo2)
# O que este comando faz?
## Fim!

ls()
#verifica os objetos criados

source( "exe1.r", echo=TRUE, print.eval=TRUE )
#mostra todos os comandos e resultados obtidos

help(source)
#mostra a pgina de ajuda com a descrio do comando
